Crossover from isotropic to directed percolation 



(N 

o 

(N 
C 

00 



43 

o 

CD 

a 



Zongzheng Zhou 1 , Ji Yang 1 , Robert M. Ziff 2 *, Youjin Deng 1 t 
1 Hefei National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei, Anhui 230027, PR China and 
2 Michigan Center for Theoretical Physics and Department of Chemical Engineering, 
University of Michigan, Ann Arbor, Michigan 48109-2136, USA 
(Dated: June 11, 2012) 

We generalize the directed percolation (DP) model by relaxing the strict directionality of DP 
such that propagation can occur in either direction but with anisotropic probabilities. We denote 
the probabilities as pi = p ■ pd and p-f = p • (1 — pd), with p representing the average occupation 
probability and Pd controlling the anisotropy. The Leath-Alexandrowicz method is used to grow a 
cluster from an active seed site. We call this model with two main growth directions biased directed 
percolation (BDP). Standard isotropic percolation (IP) and DP are the two limiting cases of the 
BDP model, corresponding to pd — 1/2 and pa = 0, 1 respectively. In this work, besides IP and DP, 
we also consider the 1/2 < pd < 1 region. Extensive Monte Carlo simulations are carried out on 
the square and the simple-cubic lattices, and the numerical data are analyzed by finite-size scaling. 
We locate the percolation thresholds of the BDP model for pd = 0.6 and 0.8, and determine various 
critical exponents. These exponents are found to be consistent with those for standard DP. We 
also determine the renormalization exponent associated with the asymmetric perturbation due to 
Pd — 1/2 near IP, and confirm that such an asymmetric scaling field is relevant at IP. 
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I. INTRODUCTION 



Directed percolation (DP), introduced in 1957 by 
Broadbent and Hammersley [l[ , is a fundamental model 
in non-equilibrium statistical mechanics and represents 
the most common dynamic universality class Q • DP has 
a very wide application, including flow in a porous rock in 
a gravitational field, forest fires, epidemic spreading, and 
surface chemical reactions Q . The DP process can be il- 
lustrated in the simple example of bond DP on the square 
lattice. Along the horizontal (vertical) edges of the lat- 
tice, the propagation occurs in a particular direction only, 
e.g., toward the right (the up). Frequently, the preferred 
spreading direction is termed "temporal," and the per- 
pendicular one is called "spatial;" the two-dimensional 
DP is thus often called "(l+l)-dimensional DP." The 
DP process has two distinct phases: the inactive phase 
for small occupation probability p where the propagation 
quickly dies out, and the active phase for large p < 1. 
Between these two phases, a transition occurs at p c . As 
the threshold p c is approached, the temporal (||) and the 
spatial (_L) correlation lengths diverge but with distinct 
critical exponents: £|| ~ \p — pd - " 11 and ^ ~ \p — p c \^" ± . 
The anisotropy is characterized by the so-called dynamic 
exponent z = V\\jvs_. For p > p Cl the order parameter 
■Poo, defined as the probability that a randomly selected 
site can generate an infinite cluster, becomes non-zero 
and its behavior can be described as Too ~ (p—Pc) , with 
/3 another critical exponent. Below the upper critical di- 
mensionality (d c + 1) with d c — 4, the three independent 
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critical exponents, v», /3, and z, are sufficient to describe 
the DP universality class. While analytical results are 
scarce for DP, even in (1 + 1) dimensions, approximation 
techniques like series expansion (341] and Monte Carlo 
simulations (9MT2I have produced fruitful results. More- 
over, after a great deal of efforts, experimental realization 
of the DP process has been achieved [l3l - [T5j in nematic 
liquid crystals, where the DP transition occurs between 
two turbulent states. 

Analogously, standard isotropic percolation (IP) [l6| is 
a fundamental model in equilibrium statistical mechan- 
ics. IP has attracted extensive research attention both 
in the physical and the mathematical communities, and 
the critical behavior is now well understood. Due to the 
isotropy, there exists only one spatial correlation length, 
which scales as £ ~ \p — p c \~ v near p c . Numerous ex- 
act results are now available in two dimensions (2D). For 
bond IP on the square lattice, the self-duality yields the 
threshold p c = 1/2 [TtJ ; the values of p c are also exactly 
known for bond and site percolation on several other lat- 
tices [2t2(|: or have been determined to a high preci- 
sion [21| . Thanks to conformal field theory and Coulomb 
gas theory [22hj25l | , the critical exponents v and fj are also 
exactly known as v = 4/3 and /3 = 5/36. 

In this work, we introduce a generalized percolation 
propagation process that contains DP and IP as two spe- 
cial cases. On a given lattice, each edge is assigned to one 
of the three possible states: occupied by a directed bond 
along a particular direction, occupied by a directed bond 
against the particular direction, or unoccupied. This is 
illustrated in Fig. Q] (a) , and the associated probabilities 
are denoted as p±, p^, and 1 — p\. — Pf , respectively. As 
a result, the percolation process has two main growth 
directions. For p^ = Pf, the symmetry between the two 
opposite directions is restored, and the system reduces 
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FIG. 1: (a) State of an edge, (b) A typical cluster in the 
BDP process. The seed is at site o, and the "infected" sites are 
denoted as solid dots. Dashed lines represent vacant bonds. 



to standard bond IP. In the limiting case p± = or 1, 
propagation against or along the particular direction is 
forbidden, and one has standard DP. We call this perco- 
lation model with two main growth directions biased di- 
rected percolation (BDP). We note that the BDP model 
is described by the field-theoretic equation in [26]. 

A natural question arises: in between standard DP and 
IP, what is the nature of phase transition for BDP? For 
later convenience, we replace parameters pi and p-\ by 
two new variables as 



Pi=P-Pd 



Pf=p.(l- p d ) 



(1) 



The parameter p is the average bond-occupation prob- 
ability (irrespective of the bond direction), and pa ac- 
counts for the anisotropy. DP corresponds to pd = or 
1, while p d = 1/2 is for IP. 

In this work, extensive Monte Carlo simulations are 
carried out for BDP in two and three dimensions. A 
dimensionlcss ratio is defined to locate the percolation 
threshold. The data are analyzed by finite-size scaling, 
and the critical exponents are determined. The numeri- 
cal results suggest that the asymmetric perturbation due 
to pd — 1/2 7^ is relevant near IP, and thus that as long 
as pd 7^ 1/2, BDP is in the DP universality class. These 
results further raise the following questions, remaining to 
be explored. For IP, is the asymmetric rcnormalization 
exponent a "new" critical exponent or related in some 
way to the known ones like v and /3? Particularly, can 
this "new" exponent be exactly obtained in two dimen- 
sions? If so, what is the exact value? 

The remainder of this work is organized as follows. 
Section II introduces the BDP model, the sampled quan- 
tities, and the associated scaling behavior. Numerical 
results are presented in Sees. Ill and IV. A brief discus- 
sion is given in Sec. V. 



II. MODEL, SAMPLED QUANTITIES, 
SCALING BEHAVIOR 

A. Model 

We shall describe in details the BDP model on the 
square lattice. The generalization to higher dimensions 
is straightforward. 

As usual in the study of DP or IP, we view the BDP 
model as a stochastic growth process, and use the Leath- 
Alexandrowicz method [13, to grow the percolation 
cluster starting from a seed site. Given the square lattice 
and the seed "o" in Fig.QJb), for each of the neighboring 
edges of site o, a random number is drawn to determine 
the edge state. If and only if the edge is occupied and 
the direction originates from the seed o, the neighboring 
site is activated and belongs to the growing cluster. For 
instance, in Fig. [TJb), the four neighboring edges of site 
o are all occupied, but site c remains inactivated because 
of the "wrong" direction. After all the four neighboring 
edges have been visited, one continues the growing pro- 
cedure from the newly added sites. In other words, one 
grows the percolation cluster shell by shell (the breadth- 
first scheme). The growth of the cluster continues un- 
til the procedure dies out or the maximum distance is 
reached, which is set at the beginning of the simulation. 



B. Sampled quantities 

In the cluster-growing process, the number of activated 
sites N(s) is recorded as a function of the shell number 
s. Let us count the shell of site "o" to be the first shell, 
the configuration in Fig. (Ub) has N — 3,5,6 for s — 
2,3,4, respectively. Besides N(s), one also records the 
Euclidean distance r of each activated site to the seed 
"o" for IP and to the y axis for the anisotropic case. The 
reason for using different definitions of r is that, for the 
anisotropic case, the average center of activated sites is 
expected to move linearly along the preferred direction, 
as s increases. Accordingly, we define a revised gyration 
radius R(s) as 



R(s) 



if N(s) = 
if N(s) > 1 



(2) 



The statistical averages Af(s) = (N(s)) and TZ(s) = 
(R(s)) are then measured, as well as their statistical 
uncertainties. We also measure the survival probabil- 
ity V{s) that at least one site remains activated at the 
sth shell and the accumulated activated site number 

>M.s) = ®Z=iNW))> 

In Monte Carlo study of critical phenomena and phase 
transitions, it is found that dimensionless ratios like the 
Binder cumulant are very useful in locating the critical 
point. Therefore, we also define a dimensionless ratio 
Qm{s) =M{2s)/M{s). 
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C. Scaling behavior 

Near the percolation threshold p c , one expects the fol- 
lowing scaling behavior 

V{s,e) ~ s- Yp F(es Y *) , 

Af(s,e) ~ s Y »N(es Y ') , 

A(s,e) ~ s YA A{es Y ') , 

TZ(s,e) ~ s YR R(es Y ') , 

Q N {s,e) ~ 2 y "Q( es y ») , (3) 

where e = p — p c represents a small deviation from p c . 
Symbols Yp 7 Yn 7 Ya,Yr, and Y e denote the associated 
critical exponents, and P, N, A, R, and Q are universal 
functions. For simplicity, only one scaling field, which 
accounts for the effect due to deviation from p c , is ex- 
plicitly included in Eq. ([3]). Right at p c , as s increases, 
the survival probability V(s) decays to zero while the 
other quantities diverge, except for the ratio Q^f which 
goes to a constant. A trivial relation is Ya = Yjv + 1. 

For standard DP (pd — or 1), exponents Yp and Yjv 
are normally denoted as <5 and 77, respectively (Y/v is also 
denoted as 9 in [2!|). It can be shown that exponent 
Y e is Y c = l/f||. Further, exponent Yr relates to 8 and 
the dynamic exponent z as Yr = —8 + 1/z, where —8 
arises from the behavior V(s) ~ s . Below the upper 
critical dimensionality (d c + 1) with d c = 4, there exist 
three independent exponents, which can be chosen as un , 
/3, and z. The others can be obtained by the scaling 
relations [29| 

vx = v\\/z, 8 = p/v\\, Ti = (<k/x-2P)/v\\, (4) 

where the last one involves the spatial dimensionality d 
and is called the hyperscaling relation. In (1 + 1) di- 
mensions, these exponents have been determined to high 
precision: v\\ = 1.733 847(6), f3 = 0.276 486(8), and z = 
1.580745(10) 0. In (2 + 1) dimensions, these exponents 
are m = 1.2890(7), (3 = 0.5812(6), and z = 1.7665(4) 

us ei nam. 

For standard IP (pd = 1/2), the shell number s is fre- 
quently called "chemical distance" [32[, accounting for 
the minimum length among all the possible paths be- 
tween the seed site and the activated sites on the sth 
shell. At p c , the length s of the chemical path relates to 
the Euclidean distance r as s ~ r dmin [H,[3J], with d m ; n > 
1 denoting the shortest-path exponent. In terms of the 
Euclidean distance r, it is known that the survival proba- 
bility scales as V(r) ~ r^^l v , the accumulated site num- 
ber A(r) ~ r 1 ! 11 ', and the p c -deviating scaling behavior 
er 1 /". This immediately yields Yp = —f3/(vd m in),YA = 
l/{vd m i n ),Y N = 7/(fd min )-l, Y R = (l-fi/v)/d min , and 
Y t = l/(^d m i n ). For IP, one has the scaling relation as 

j/u = d- 20/v. (5) 

In 2D, v and (3 are exactly known as v = 4/3 and 
(3 = 5/36, which yield 7/1/ = 43/24 1.79166... and 



/3 jv = 5/48 ps 0.104166.... The shortest-path ex- 
ponent d m in, together with the so-called backbone ex- 
ponent, is among the few critical exponents of which 
the exact values are not known for the 2D percolation 
universality class. It was conjectured to be d m i n = 
217/192 = 1.13020... (35|, and some recent estimates 
are 1.1306(3) [H and 1.13078(5) H3- In three dimen- 
sions, no exact results are available, and the numerical 
estimates are d min = 1.374(4) [H, fi/v = 0.4774(1) and 
v = 0.8734(6) j3|, which yield (3 = 0.4170(4). 



III. RESULTS 

In this work, we consider the BDP model on the 
square lattice for 2D and the simple-cubic lattice for 
3D. The simulation applies the aforementioned Leath- 
Alexandrowicz growth method. The dimensionless ratio 
is used to locate the percolation threshold p c . Ac- 
cording to Eq. ([3]), ratio is expected to have an ap- 
proximate common intersection at p c for different shell 
number s. At the threshold p c , as s — } 00, the common 
intersection converges to a universal value 2 Yn and the 
slope of Qjv" increases as s Yc . 



A. Standard IP 

Standard IP corresponds to pd = 1/2. Monte Carlo 
simulation was carried out up to s max = 8192 for 2D 
and 2048 for 3D. About 10 8 samples were taken for each 
data point on each lattice. The Q^f data are shown in 
Fig. [5] Indeed, we find an approximate common intersec- 
tion near p = 1 and 0.4976 for 2D and 3D, respectively. 
This agrees with the known threshold p c /2 = 1/2 (2D) 
and 0.248 812 6(5) (3D) [H. Note that, since the occu- 
pied bond can propagate the growth process only if it 
has the correct orientation, there is a factor of 2 differ- 
ence between the bond-occupation probability p here and 
the p of the equivalent bond percolation probability. 

To have a better estimate of p Cl according to a least- 
squared criterion, the Q_\f data are fitted by 

4 

Qm(s, e) = Qat, c + qke k s kY ' + b lS yi 
fe=i 

+b 2 s-' 2 + ces Y ' +yi + ne 2 s Y ' + ... , (6) 

which is obtained by Taylor-expanding Eq. (J3j> and tak- 
ing into account finite-size corrections due to the leading 
irrelevant scaling field and analytical background contri- 
bution. These are described by the two terms with ampli- 
tudes b\ and 62 , of which the term with n arises from the 
nonlincarity of the relevant scaling field in terms of the 
deviation e, and the one with c accounts for the combined 
effect of the leading relevant and irrelevant scaling fields. 
In the fits, various formulas are tried, which correspond 
to different combinations of those terms in Eq. ([5]) . For a 
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FIG. 2: Ratio Qm for IP in 2D (top) and 3D (bottom). 



given formula, the Q^ data for small s < s min are grad- 
ually excluded from the fits to see how the residual \ 2 
changes with respect to s mm - The results from different 
formulas are compared with each other to estimate the 
possible systematic errors. In two dimensions, we obtain 
Pc = 1.000 000(4), Q u ,c = 1.499 5(1), Y e = 0.664(3), and 
j/i = —0.96(6). Note that the leading irrelevant thermal 
scaling field is u> — —2 for 2D percolation universality 



411 ] : apparently, the leading correction exponent yi = 
—0.96 does not correspond to oj. Instead, y\ should be 
associated with the chemical distance. From the relations 
Qm.c = 2 Yjv , Y N = 7/(^rf min ) - 1, and Y e = \/{vd min ), 
and the exact values j/f = 43/24 and 1/V = 3/4, we 
determine d min = 1.130 76(10) from Qx >c = 1.499 5(1), 
and d min = 1.130(6) from Y e = 0.664(3). ' 

In three dimensions, our results are p c = 0.497 624(1), 
Q M c = 1.400(1), Y e = 0.830(1), and y x = -0.8(2). Our 
estimate of p c /2 = 0.248 812 0(5) agrees with the exist- 
ing one 0.248 812 6(5) [40l| . and has a comparable error 
margin. 

To estimate other critical exponents, we simulate right 
at the threshold p/2 = 1/2 for 2D and 0.248 812 
for 3D. The simulation was carried out for s up to 
s max = 8192 for 2D and 2048 for 3D. Further, to elimi- 
nate one more unknown parameter in the fits, we mea- 
sure the dimensionlcss ratios Qv(s) — P(2s)/V(s) and 
Qti = lZ(2s)/1Z(s). These Q data are fitted by 



Q(s) = Q c + b lS yi + b 2 s 



(7) 



In two dimensions, the results are Qv,c = 0.9382(1) and 
yi = -0.80(7) for Q V: and Q n .c = L7318(2) and yi = 
—0.9(1) for Qti. For all these three ratios, the leading 
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FIG. 3: Ratio Qm for standard DP in 2D (top) and 3D (bot- 
tom). 



correction is described by an exponent y\ « —1. Taking 
into account the exact values (3 /v = 5/48, one has <i m i n = 
1.132(2) from Qv,c and d m - m = 1.130 7(3) from Qn,c- 

In three dimensions, the results are Qj>,c = 0.7865(2) 
and yi = —0.7(2) for Q-p, and Qn,c = 1.3020(3) and yi = 
—0.9(2) for Qti- Combining the estimate Qv,c and Qn,c 
together, one has (3/v = 0.4765(8) and d m in = 1-375(1). 
Our result rf min = 1.375(1) agrees well with the existing 
result d m i n = 1.374(4) [38(, and significantly improves 
the error margin. 

For comparison, these results are summarized in Ta- 
ble D 



B. Standard DP 

We simulate standard DP by taking pd = 1. The sim- 
ulation was carried out for s up to s max = 16384 for 2D, 
and 2048 for 3D. The number of samples for each data 
point is about 8 x 10 8 in 2D and 1.6 x 10 8 in 3D. 

The Qjsf data are shown in Fig. [3] A good intersec- 
tion is observed for both 2D and 3D, which yields p c ~ 
0.64470 for 2D and 0.38222 for 3D, from a rough visual 
fitting. We fit the Qn data more precisely using Eq. © . 
On the square lattice, we obtain p c = 0.644 700 5(8), 
Qm,c = 1-242 9(2), Y t = 0.576(3), and y Y = -0.9(1). The 
estimate of the percolation threshold agrees well with the 
existing more precise result 0.644 700185(5) 0- From 
the relations Qa/\c = 2 Yn = 2 V and Y t = we have 

r/ = 0.313 7(2) and z/y = 1.736(9). On the simple-cubic 
lattice, our results are p c = 0.382 225 6(5), Y t = 0.777(2), 
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^min 


Pc/2 


2D (known) 
(present) 


5/36 [16, 22-25] 
0.138 7(10) 


4/3 [16, 22-25] 
1.332(6) 


1.130 6(3) [35-37] 
1.130 76(10) 


1/2 [16, 17J 
0.500 000(2) 


3D (known) 
(present) 


0.4167(4) 
0.417(1) 


0.873 4(6) [39] 
0.876(2) 


1.374(4) [38] 
1.375(1) 


0.248 812 6(5) [40] 
0.248 812 0(5) 



TABLE I: Percolation thresholds and critical exponents for IP. 
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FIG. 4: Ratio Qm for BDP in 2D. The top (bottom) panel 
corresponds to pd = 0.8 (0.6) case. 



FIG. 5: Ratio Qm for BDP in 3D. The top (bottom) panel 
corresponds to pd = 0.8 (0.6) case. 



C. BDP 



Qm.c = 1.1738(1), which yield 77 = 0.2312(1) and v\\ = 
1.287(4). Here the y\ is too small to estimate since the 
numerical data of s > 24 can be well described even 
though we do not include any corrections. The agreement 
of p c with the existing estimate p c = 0.382 224 64(4) [U 
is within two standard deviations. 

Analogously, we simulate right at the percolation 
threshold p c = 0.644 700 185 for 2D and p c = 0.382 224 64 
for 3D. The dimensionless ratios Q-p and Qn are mea- 
sured, and the data are fitted by Eq. ([7]). For 2D, 
the results are Q v , c = 0.89537(5), y x = -0.98(5) and 
Qn,c = 1.3882(1),' y x = -1.1(1), which yield Y P = 
8 = 0.15944(9) and Y R = (-S + 1/z) = 0.47322(10). 
Taking into account the estimates of z/11 and 6, one has 
v±_ = 1.098(6). For 3D, the results arc Q P , C = 0.7311(4) 
and Qn,c = 1.0822(1), which yield that 5 = 0.4519(8) 
and V± _ = 0.728(4). 

These results are listed in Table ILTl 



For the purpose of studying BDP, we choose pd = 0.6 
and 0.8. The simulation was carried out for s up to 
s max = 16384 for 2D and 2048 for 3D. About 2 x 10 8 
samples were taken for each data point in each case. 

The Q_\f data are shown in Fig. U] for 2D and Fig. [5] 
for 3D. The transitions are also clearly observed, but the 
approximate common intersections are not as good as 
those for standard DP and IP. This suggests the existence 
of additional finite-size corrections. 

The data are also fitted by Eq. © according to a 
least-squared criterion. To account for the possible exis- 
tence of additional corrections, we replace the terms in 
Eq. ©, with 61, 6 2 , and c, by b, s ! «+M yi + ces« i+v *. The 
exponent y\ is fixed at —1, in accordance with our above 
estimate of y\ for both standard IP and DP. Indeed, the 
new source of finite-size correction can be identified in 
the fits, which yield yi = —0.5(2) both in 2D and 3D. 
The results for p c , rj = log 2 Qm,ci and v\\ = 1/Y e are 
summarized in Tabic HT1 

The determination of the critical exponents S and z 
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Ref. 
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TABLE II: Percolation thresholds and critical exponents for standard DP (pd = 1) and BDP (pd < 1). The numbers in the row 
with reference are the existing results. Clearly, standard DP and BDP with pd = 0.8, 0.6 share the same critical exponents. 
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FIG. 7: Top: Ratio Qm versus (p d - 1/2) 2 with p — 1 on 
square lattice. Bottom: Log-log plot of pdc — 1/2 versus 1 — p c 
for the transition line (pdcPc) near IP. The dashed line has 
slope 0.754. 



is obtained in an analogous way by simulating at the 
estimated percolation threshold, and the results are listed 
in Table HH 

The results in Table ITTI strongly suggest that, as long as 
Pd deviates from 1/2, the system falls into the standard 
DP universality class. For an illustration, we make the 
log-log plot of the critical quantity TV" versus the shell 
number s in Fig. |51 Clearly, the slope for pd = 1/2 is 
distinct from those for the other cases, which are inde- 
pendent of Pd- 

IV. CROSSOVER EXPONENT 

The fact that BDP for pd ^ 1/2 is in the DP universal- 
ity means that in the language of renormalization group 



theory, the operator associated with the asymmetric per- 
turbation is relevant near the IP fixed point. To confirm 
this, we simulate BDP near IP with p = p c = 1 by vary- 
ing ed — Pd — 1/2. The simulation is up to s max = 8192, 
and ed is set at 0, 10 -3 and 2 x 10 -3 . The results for 
Q_\f in two dimensions are shown in Fig. [7] versus e 2 d ] 
note that BDPs for ±td are identical. These Qj\f data 
are also analyzed by Eq. ^ with Y e being replaced by 
the exponent Y £d for the symmetric scaling field and the 
odd terms with respect to ed being set zero. We obtain 
Y Cd = 0.500(5), which suggests that Y ed may be exactly 
1/2. 

According to scaling theory, the phase transition line 
(PcPdc) approaches to the critical IP (p c = l,pdc = 1/2) 
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FIG. 8: Phase diagram of the BDP model in 2D (left) and 3D 
(right). The pd — 1/2 line corresponds to isotropic percola- 
tion. The diagram for pd < 1/2 is drawn by symmetry. The 
arrows represent the direction of the renormalization flows. 



as m 



l-p c oc \(p dc -l/2)\^ 



(8) 



where <j) — Y Cd /Y c is the so-called crossover exponent. 
We carried out some Monte Carlo simulations and deter- 
mined a set of critical points near IP; they are 13 critical 
points with p dc £ [0.52,0.6]. In Fig. [71 we plot p dc — 1/2 
versus 1 — p c in log-log scale, which indeed has slope ap- 
proximately equal to <j) = Y ed /Y e = 0.754. 

We also perform a similar study near the critical IP in 
3D, and obtain Y td = 0.56(1) and = 0.67(1). 



V. DISCUSSION 

We introduce a biased directed percolation model, 
which includes standard isotropic and directed percola- 
tion as two special cases. Large-scale Monte Carlo simu- 
lations are carried out in two and three dimensions. We 
find that the operator associated with the anisotropy is 
relevant near the IP fixed point, which implies that BDP 
in the region p d ^ 1/2 is in the DP university class. On 
this basis, the phase diagram and the associated renor- 
malization flows are shown in Fig. [8] Since the upper crit- 
ical dimensionality is different for standard IP and DP, 
it is not clear whether the similar renormalization flows 
would hold in higher dimensions. We mention that such 
crossover phenomena have attract much attention both 
in the fields of equilibrium and non-equilibrium statisti- 
cal mechanics [43j-|49(. In retrospect, it is not surprising 
that the asymmetric perturbation is relevant near IP. At 
IP, all the directions are equivalent and "spatial" and 
"temporal" directions cannot be defined. However, as 
soon as p d — 1/2 ^ 0, such a symmetry is broken and 
the center of the activated sites moves along the "tempo- 
ral" direction as the growing process continues. It is also 



plausible that as long as the "spatial" and "temporal" 
symmetry is not restored, such an asymmetric pertur- 
bation is irrelevant near DP. This is similar to the fact 
that asymmetric diffusion on the basic contact process 
is irrelevant [49j . In terms of the chemical distance s, 
the effect from the anisotropy can be asymptotically de- 
scribed as oc (p d - l/2)s Y *i with Y €d (2D) = 0.500(5) 
and Y £d (3D) = 0.589(10). One can also use the Eu- 
clidean distance r to describe such an anisotropic effect 
as oc {p d — \j2)r 1 l Vd with Y td = 1 / '(i/^b) • Substituting 
the d mm value into Y ed , one obtains v d {2D) = 1.77(1) 
and v d (ZD) = 1.30(2). 



When viewing standard isotropic percolation in the 
framework of BDP, one observes that two independent 
critical exponents, e.g., v and /3, are no longer suffi- 
cient to describe the critical scaling behavior. In this 
case, the shortest-path exponent e? m i n appears naturally 
and becomes indispensable, and thus isotropic percola- 
tion also has three independent critical exponents. Our 
estimate of d m in significantly improves over the existing 
results both in two and three dimensions. Our result 
rfmin = 1-130 76(10) does not agree with the recently 
conjectured value 217/192 (3f| in two dimensions. This 
result appears to refute the conjectured value. On the 
other hand, we note that, in terms of the chemical dis- 
tance s, a new source of finite-size corrections occurs in 
the scaling behavior, and these corrections are not well 
understood yet. Further, we observe that the restored 
symmetry for IP can be regarded as v\\ = in the BDP 
model. In some cases, the coincidence of two critical 
exponents may suggest the existence of logarithmic cor- 
rections of the log or log log form, and they can be either 
additive or multiplicative. In practice, logarithmic finite- 
size corrections have indeed been observed for standard 
isotropic percolation in two dimensions 2l|, which is in 
terms of Euclidean distance. In this sense, we cannot 
entirely exclude the possibility that the tiny difference 
between the present numerical result for d m ; n and the 
conjectured value arises from some unknown corrections 
that have not been taken into account in the numerical 
analysis. Numerical investigation of this problem seems 
very difficult if not impossible. Nevertheless, since the 
exact value of <i m i n is conjectured as a function of q for 
the g-state Potts model |35| . one can accumulate more 
numerical evidence by studying the q ^ 1 case. 



Finally, the numerical estimate of the critical exponent 
Y ed or v d due to the asymmetric perturbation near IP 
raises a question: is it a "new" independent critical ex- 
ponent or simply related in some way to the known ones 
like f3, v, and d m i n ? In particular, in two dimensions, one 
would ask whether v d or Y £d can be exactly obtained in 
the framework of Stochastic Loewner Evolution (SLE), 
conformal field theory or Coulomb gas theory. 
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